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(— I . Abstract 



The nonlinear system of equations of relativistic Lorentz electrodynamics (LED) is stud- 
ied in a "gyroscopic setup" in which the Lorentz electron is assumed to remain at rest, 
. leaving the electromagnetic fields and the particle spin as the only dynamical degrees of 

^ I freedom. The global existence and uniqueness of this gyroscopic spin-plus-field dynamics 

^ ' in unbounded space is proven. It is further shown that for rotation-reflection symmetric 

initial data any gyroscopic solution also satisfies the world-line equations consistent with a 
non-moving Lorentz electron, thus furnishing a proper solution of the complete set of equa- 
tions of LED. Rotation-reflection symmetric scattering is shown to occur for sufficiently 
small ratio of electrostatic to (positive) bare rest mass, with deviations from the station- 
ary spin state dying out exponentially fast through radiation damping. The previously 
proven result that the renormalized spinning Lorentz electron evolves like a soliton in scat- 
tering processes combined with the present results that scattering does occur establish the 
solitonic character of the renormalized Lorentz electron. 

To appear in : Letters in Mathematical Physics 
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1 Introduction 



In recent years the century-old Lorentz program of electrodynamics |jTO| has attracted the at- 



tention of mathematical physicists. Most of the rigorous results established so far belong to the 



semi-relativistic Abraham model and are surveyed in [|1J]. Recently [|1| the authors presented 
the first properly renormalized approach to truly relativistic Lorentz electrodynamics (LED), 
picking up on the pioneering work . While the dynamical initial value problem for the model 
in 1^ is seriously singular, our Lorentz-covariant LED displays most of the features considered 
crucial for a realistic, consistent classical electrodynamics, namely: 

• the Cauchy problem for the evolution of the physical state in "massive" LED with strictly 
positive bare rest mass and bare moment of inertia is regular; 

• the pre- and post-scattering values of the renormalized electron rest mass and electron spin 
magnitude are identical, i.e. the Lorentz electron evolves as soliton in scattering processes; 

• there exists a simple curve in the charge - bare rest mass - radius - gyration frequency — 
parameter space of the stationary bare particle on which the stationary renormalized particle 
data charge, magnetic moment, and mass match the empirical electron data without involving 
superluminal gyration speeds. 

In Ijl]] we also studied LED's renormalization flow to vanishing bare rest mass with empirically 
matched data when the positive bare mass and charge are distributed on the surface of a sphere. 
The renormalized "purely electromagnetic" LED which emerges in the limit has the following 
additional characteristics: 

• the renormalized purely electromagnetic LED constitutes a classical field theory equipped 
with an ultraviolet cutoff at about the physical electron's Compton length; 

• in the limit of vanishing bare rest mass the equatorial gyration speed reaches the speed of 
light and the bare gyrational mass converges to a "photonic" mass; 

• in the same limit, the renormalized spin magnitude converges to 3^/2, up to corrections of 
order a (Sommerfeld's fine structure constant). 

In this paper we supply several rigorous results regarding scattering in LED conditioned on 
a straight particle world-line. For a straight particle world-line the set of Maxwell-Lorentz plus 
gyrational equations decouples from the world-line equations, which in turn become constraint 
equations that impose on the admissible set of initial conditions for the dynamical equations for 
the spin and the electromagnetic fields. We prove that all physically reasonable Cauchy data 
for the fields and the spin launch unique global forward and backward evolutions of Maxwell- 
Lorentz plus gyrational equations. For rotation-reflection symmetric data it is shown that these 
gyroscopic solutions satisfy the world-line equations for a non-moving Lorentz electron, too, 
thus furnishing solutions of the complete set of equations of LED. Rotation-reflection symmetric 
scattering is shown to occur if the ratio of electrostatic to bare rest mass is smaller than 1. 
The previously proven result that the renormalized spinning Lorentz electron evolves like a 
soliton in scattering processes combined with the present results that scattering does occur 
establish the solitonic character of the renormalized Lorentz electron. It is further shown that 
the rotation-reflection symmetric deviations from the soliton state die out exponentially fast 
through radiation damping. The results proven here are somewhat stronger and cover more 
general mass and charge densities than announced in 
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2 Notation 



We use the notation of which largely follows the conventions of 0. Abstract Minkowski space 
is identified with M^'^, equipped with a Lorentzian metric of signature +2. Thus, any orthonormal 
basis {eo, ei, 62, 63} satisfies the elementary inner product rules Cq ■ eo = —1, ■ = 1 for 
/X > 0, and e^^ ■ = for fi ^ u. A constant basis defines a Lorentz frame, denoted JFl. We 
use X, y, etc. to denote four-vectors representing events in spacetime. With respect to JF^, 
we decompose x into time-plus-space components thus, x = {ct,x), where x = {x^,x'^,x^) is 
a "point in space," and t = x^/c an "instant of time," where c is the speed of light in vacuo. 
Henceforth we shall use units in which c = 1. We call v spacelike, lightlike, or timelike according 
as V • V > 0, V ■ V = 0, or V ■ V < 0, respectively. We define ||v|| as the principal value of 
(v ■ v)^/^. The tensor product is defined by its inner-product action on four- vectors 
thus, (e^ By) ■ c = e^{ey ■ c) and c ■ (e^ ® e^,) = (e^ ■ c)ey. In general a rank-two tensor reads 
T = T^'^e^ ® e^, and if T^'^ = ±T^^ it is symmetric (+ sign), respectively anti- symmetric (— 
sign). The metric tensor g = g^^e^ ® e^,, with g^^ = ■ e^, is clearly symmetric and has the 
same components g'^" in all Lorentz frames. Notice that g acts as identity on four-vectors, i.e. 
g ■ V = V. A particular class of anti-symmetric tensors is given by the exterior product between 
two four- vectors, aAb = a^b—b^a. Finally, [A, B]_|_ = A ■ B ± B ■ A is the (anti-) commutator 
of any two tensors of rank two A and B. 

For a differentiable function /(x) we denote by Vg/ its four-gradient w.r.t. g. In time- 
plus-space decomposition, Vg/(x) = (— (9^^/, V/), where V is the usual three-gradient. The 
four-curl of a differentiable four-vector function is defined in analogy with the conventional curl 
as the anti-symmetric four tensor function 

VgAA(x) =£^^^V^®e,(eA- Vg)(e,- A) (2.1) 

where the the entries of the rank-four Levi-Civita tensor. The four-Laplacian with 

respect to g is just the (negative) d'Alembertian, or wave operator, i.e. Ag = Vg ■ Vg = — □. 

3 Covariant massive LED with a straight particle world-Une 

In this section we present the manifestly covariant equations of massive LED [|T| for the special 
case that the particle's world-line is straight. 

3.1 Kinematical pre-requisites 

We recall that the particle's world-line is a map r x = q(r), where dr = V— dx ■ dx, with dx 
taken along the world-line, is the invariant proper-time element. The map r u = q(r), where 
q is the particle's four-velocity, is the world-hodograph. The world- gyrograph of the particle is 
an anti-symmetric tensor- valued map r ^e{t) of space-space type with respect to u (i.e. 
Oe ■ u = 0) which describes the angular velocity of the inert gyrational motions of the Lorentz 
particle that may occur in excess of the inertia-free Thomas precession. Thus, = fl — u A u, 
where fl is the angular velocity tensor of the particle's co-rotating body frame, while u A u is 
the familiar angular velocity tensor of Fermi- Walker transport 0. 

For a straight world-line ^(t) = uot + Qq the particle's four- velocity is a constant four- vector, 
u(r) = uq for all r. A constant four- velocity in turn implies that fig = f2. 
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3.2 Field equations 



The electromagnetic Maxwell-Lorentz fields are gathered into the anti-symmetric rank-two Fara- 
day tensor field x F(x), which satisfies the manifestly covariant Maxwell-Lorentz equations 

Vg-'^F = 0, (3.1) 



Vg-F = 47rJ, (3.2) 

where *F is the (left) Hodge dual of F and J is the charge- current density four- vector field, given 
by Nodvik's M manifestly covariant expression 



/ + 00 
(uo - f^(r) ■ x) /e( ||x - q(r)|| ) 6{u, ■ (x - q,) + r) dr , 
-oo 



(3.3) 



where fe '■ [0,R] — ^ is the SO (3) invariant charge "density" of the Lorentz particle, and 
< i? < cx) its radius. For a Lorentz electron, f^^ /e(|a;|)d^x = — e, where e > the elementary 
charge. Conditioned on the world-line r i— >• q(r) = uqt + Qq and gyrograph r h-* ^e{t) = ^{'t) 
being given, the Maxwell-Lorentz equations are linear equations for F. 

3.3 World-gyrograph equations 

The equations for the gyrograph are 

^S, = t, (3.4) 



where 



Sb(r)=/ {y-uoT)A^^M=^=f^{\\y-uoT\\)6{uo-y + T)d'y (3.5) 



is the anti-symmetric tensor of bare Minkowski spin (about q(r) = uqt + q^) associated with 
the gyrational motion of the S0{3) invariant bare rest mass "density" /m : [0,-R] — ffi"*" of the 
particle, while 

t(r) = [ (y - uor) A (F(y) • (uo - f^(r) ■ y))^/e( \\y - uot\\ ) 6{uo ■ y + r)d^y (3.6) 

is the Abraham-Lorentz type Minkowski torque, with = (g + Uq ® uo) ■ ^■ 

3.4 World-line equations 

The world-line equations are 

Ap = f, (3.7) 

where 

p(r) = M(r) • Uo (3.8) 
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is the Minkowski momentum four-vector of the particle, with M = + Mbg its symmetric 
Minkowski tensor mass, where 

MnW = - / [{y - uqt) ® (y - uqt), [F(y), ri^M] J /e( \\y - uo^H ) 5(uo • y + T)d^y 

(3.9) 

is the Nodvik tensor mass extracted from the Minkowski momentum four- vector associated 
with electromagnetic spin-orbit coupling given in ||^, and where 

Mb(r)= / (l-\\n,-yfy^U{\\y-Uor\\)6{uo-y + r)dS (3.10) 
is the gyrational bare mass Finally, 

f{r)= [ Hy)-{uo-^.{r)-y)my-UoT\\)6{uo-y + T)d'y (3.11) 
is the Abraham-Lorentz type Minkowski force . 

4 The Cauchy problem for the state in LED 

We now choose a convenient Lorentz frame, called the "laboratory frame" J-'iah, in which the 
space-plus-time decomposition of our manifestly covariant equations takes a simple form. In 
particular, since we consider only evolutions for which u(r) = uq for all r, we can work with 
the standard foliation of space-time in our frame J-'\ah- The standard foliation of J-'i^h consists 
of the level sets Tjr^^^^x) = t of the function T;Fj^^(x) = — eo ■ x, which has a constant timelike 
four-gradient VgT:ri^b(^) — ~^o- The space-plus-time decomposition of events in jFjab written 
as {t,x), is understood w.r.t. this standard foliation. 

By a boost we can achieve that the timelike unit vector eg = (1,0,0,0) of J-'i^b coincides 
with the four-velocity of the particle, i.e. uq = eo. By at most a spacetime translation we 
can furthermore assume that q(0) = in J-'iah, so that the particle's space position is at the 
origin of the space hypersurface of J-'\s.h, and that laboratory time t and particle proper-time r 
coincide. Accordingly, from now on we will write t in place of r. The world-line as seen in J-iab 
is now simply given by q{t) = (t,0). As for the gyrograph, since = we will henceforth 
simply omit the subscript e- In J-'iah we clearly have f2(t) ■ eo = for all t, so that f2 is dual 
to a spacelike four-vector w(t) which satisfies fl{t) ■ w{t) = and w{t) ■ eg = for all t. 
Hence, in our jFi^b we have w = {0,uj), where a;(t) is the usual angular velocity three- vector, 
directed along the instantaneous (i.e., at time t) axis of body gyration in the space hypersurface 
of J-\a.h- Finally, the field tensor F(x) at x is decomposed as usual into its electric and magnetic 
Maxwell-Lorentz components w.r.t. the standard foliation of J-'iah, here conveniently grouped 
together as a complex electromagnetic three- vector field, 

gix,t)=£{x,t)+tBix,t), (4.1) 

whose real and imaginary part are, respectively, the electric (i.e. time-space) and magnetic (i.e. 
space-space) components of the field tensor F in J-'iah- Since by hypothesis q(t) = (t, 0) for all t, 
the state at time t in LED is uniquely characterized by specifying uj(t) and G{- ,t). 
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4.1 Evolution equations 

The covariant equations now decompose into a system of first-order evolution equations for the 
state variables of LED, plus a set of constraint equations. We begin with the evolution equations. 

4.1.1 Field equations 

Beginning with the covariant field equations, we note the space-plus-time decomposition of the 
current density four-vector as J(x) = (1, x x)fe {\x\). The space components of the covariant 
field equations combine into the Maxwell-Lorentz evolution equations for Q, 

dtG{x,t) = -iV X g{x,t)-ATiuj{t) xxU{\x\) , (4.2) 

where dt means first-order partial derivative w.r.t. Lorentz time and Vx is the standard curl 
operator. 

4.- 1-2 Spin equations 

Turning next to the gyrational equations, we recall that fig is dual to the space vector uj. In the 
same vein, the space projector g+ uq CS) uq under the integral in (|3^ ) guarantees the space-space 
character of Sb w.r.t. uq, i.e. Sb ■ uq = for all r, so that the bare spin Minkowski tensor ( |3.5D 
is dual to the space vector of bare spin, 

f XX (i^(t) X a;) „ o , , 

Sb(t)=/ -^=^=^=I=J=U\x\)d'x, (4.3) 



^K" — \0j{t) X X\'^ 

and the Minkowski torque ( p.6| ) is dual to the torque space vector 

t{t)= XX (£{x,t) + {u;{t)xx)xB{x,t))fe{\x\)d^x. (4.4) 

Equation (|3.4|) together with (|3.6|) is therefore dual to the evolution equation 

^Sb = *, (4.5) 
for uj{t). 

4.2 Constraint equations 

4-2.1 Divergence equations 

The time components of the covariant field equations combine into the Maxwell-Lorentz diver- 
gence equation 

V-g{x,t)=A7Tf,{\x\). (4.6) 



Notice that ( |4.6| ) is merely a constraint on the set of initial data, for the (three-) divergence 
of ([4.2|) implies that a solution Q{x,t) of ([4.2|) for given uj{t) x a;/e(|a;|) automatically satisfies 
( [4.6|) for alH > if the initial data Qq = Sq + iBo satisfy the constraint ( [4.6|) at time t = 0, i.e. 
if V-C?o(a^) =47r/e(|a^|). 



6 



4-2.2 World-line equations 

The four-momentum p has the space-plus-time decomposition p = (Mb, Ng ■ cj), where 



Mb(t)=/ „ :/m(|a;|)d^x (4.7) 

1 CL?(rj X x\ 



is the bare gyrational mass at time t, and where 

Ne(t) = / a;® fa; x^(a;,t))/e(|a;|)d3x (4.8) 

is a spin-orbit couphng tensor. Furthermore, the Abraham-Lorentz type Minkowski force now 
has the space-plus-time decomposition f = (P, /), where 

P{t) = uj{t)- f (xxE{x,t))U\x\)d?x (4.9) 

is the power delivered by the field to the particle, and where 

f{t)= f (£{x,t) + {u{t)yix)yiB{x,t))U\x\)d^x (4.10) 

is the Abraham-Lorentz force on the particle. The space-plus-time decomposition of the world- 
line equation then becomes 

1a4 = P (4.11) 



and 



^(Ne-^) = /. (4.12) 



Despite their appearance, equations ( [4.11| ) and ( [4.12| ) are not evolution equations for the world- 



line; instead, they have to be satisfied by the active state variables a;(t) and Gi-,t) to ensure 
consistency with the constraint that the world-line is given by q{t) = e^t in J^i^h- However, 
we shall show that ( [4.11 ) is automatically satisfied for all time by any solution of the evolution 



equations for spin and fields that obeys the divergence equations initially. This leaves ( |4.12| ) as 



the only true constraint equation coming from the world-line equation. While we will show that 
certain symmetric initial conditions launch a dynamics consistent with ( |4.12| ), it seems difficult 



to precisely characterize the complete set of initial conditions that will launch such a consistent 
dynamics. 

4.3 Cauchy data 

The field evolution equation ( |4.2| ) are supplemented by initial data consistent with the constraint 
equations (|4.6| ) and satisfying the asymptotic condition that Q{x,t) ^ as \x\ —>■ oo, the real 
part as S{x,t) ~ — ea:;/|a:;|^ + o(|a;|~^), the imaginary part satisfying \B\ = 0{\x\~^). 

Equation (|4.5| ) is to be supplemented by initial data uj{0) = ujq satisfying the requirement of 
strict subluminality, |<-Jo|-R < 1, or subluminality, \ujq\R ^ 1, depending on the choice of f^. 
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Viewed from a dynamical systems perspective, Cauchy data may be prescribed in any consis- 
tent manner, and for our existence and uniqueness result of a strong solution in some weighted 
V" norm we only need that the cumulative time integral of the wave fields over the support of 
the particle stays bounded. However, the scope of LED as a theory, in the classical limit, of the 
dynamics of an electron coupled to the electromagnetic fields, its self-fields included, basically 
limits the physically sensible choices of initial data to a stationary electron well-separated from 
some localized radiation field that has compact support in space disjoint from the fixed support 
of the particle. To have a dynamically interesting scenario, the time-evolved support of the 
initial radiation fields should eventually overlap with the support of the electron. 



5 Gyroscopic LED 

We study first the subsystem of equations obtained by neglecting the world- line equations ( [4.11| ), 
( [4.12| ) from the LED with a straight world-line. For obvious reasons, we will call this dynamical 
model the gyroscopic LED. 

We have to solve the Maxwell-Lorentz equations ( [4. 21 ), ( [4. 61 ) for the field ( |4.1| ) together with 



the gyrograph equations ( ^4.5| ), ( [4. 4]) for the bare spin ( |4.3| ). Our strategy is to solve first the 
Maxwell-Lorentz equations in terms of integral representations involving the unknown bare spin 
dynamics. Inserting this representation into the gyrograph equation, we rewrite the latter into a 
fixed point problem for Sb(t). We then prove that the fixed point map is a Lipschitz map, from 
which the global well-posedness of the gyroscopic problem follows. Subsequently we will show 
that the gyroscopic problem conserves the energy, angular momentum and the canonical spin 
magnitude, but generally not the linear momentum. Energy conservation is coincidental with 



the fact that (|4.11|) is automatically satisfied by a gyroscopic solution. 



5.1 Forward integration of the Maxwell-Lorentz equations 

We recall that in virtue of the homogeneous Maxwell-Lorentz equations ( p.l|) , there exists a 
(non- unique) four- vector field A satisfying the Lorentz gauge Vg ■ A = such that F = Vg A A. 
The inhomogeneous Maxwell-Lorentz equation ( |3.2| ) then becomes the inhomogeneous wave 
equation □A(x) = 47rJ(x). Recalling furthermore the time-plus-space decomposition for the 
current density four-vector, J(x) = (l,a;(t) x a;)/e(|a;|), and introducing the time-plus-space 
decomposition for the electromagnetic potential four-vector as A(x) = (0(a;, t),.4(a;, t)), the 
equation F = Vg A A becomes 

g{x, t) = -V(f){x, t) - dtA{x, t)+iV X Aix, t). (5.1) 

The Coulomb potential (p and vector potential A. satisfy the inhomogeneous wave equations 

Dcf>{x,t)=47if,{\x\), (5.2) 



DA{x,t) = A7if^{\x\)uj{t) X X, (5.3) 

supplemented (i) by the asymptotic conditions (f){x,t) ~ — e|a3|~-'^ and A{x,t) HqX x \x\~^ as 
— i> CO, for all t G M, where /Xq is the particle's magnetic moment at t = 0, 

f^o = l [ XX {ujqX a;)/c(|a;|)d^a; (5.4) 
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with ujq = tt'(O), and (ii) by compatible Cauchy data at t = 0. 

We first integrate the wave equations (^.2|) , { f).3\) for the potentials and A. Clearly, (^]2| 
is solved by 4>{x,t) = (f)coni{x) + 0wave(ic,t), where 







/r3 \x - y\ 





t)=[ Lo + e(t -\x- y\){uj{t -\x- y\) - cuq)) x ^^/c( \y\ )d^y (5.6) 
Jr3\ 'J \x-y\ ^ ' 



Coui(a;)= / . rfJ\y\)d'y (5.5) 



is the static Coulomb potential for /e and 0wave(a^,^) is a solution of the homogeneous scalar 
wave equation □(/)wave(3^, t) = 0. After at most a gauge transformation, we may assume that 
0„ave = 0. Next, ( |0| ) for t > is solved by A{x, t) = ^source(2;, t) + Awa.ve{x, t), where 

•'^-source (''^; 

solves the inhomogeneous vector wave equation (|5.3| ) (O is the Heaviside function), and where 
Awa.ve{x,t) solves the homogeneous vector wave equation DAwa.ve{x,t) = for initial data 
>iwave(a?,0) = Aq{x) and dtA„^^e{x,0) = -£'q{x), where ^o(^) = -^oix) - Aource(a;, 0), with 
^o(^) the initial magnetic vector potential, and where Sq{x) = Sq{x) + V0coui(^c), with Sq{x) 
the initial electric field strength. Thus ^wave is given by Kirchhoff 's formula 

Xave(a;,t) = -7 / ^'Md% + wXl I -^M^^) , (5-7) 

JdBt(x) Ot\t JdBt{x) ' 

where df2y is the uniform surface measure on dBtix) divided by 47r. 

5.2 Canonical form of the gyrograph equation 

With the help of the potential representation of Q we now rewrite ( [4. 51) , (|4.4|) into the more 
accessible canonical format. Recalling that £{x,t) = — V0coui(^c) — dtA{x,t), with (j)cou\{x) 
given in ( |5.5|) , and with B{x,t) = V x A{x,t), and noticing that x x V0coui(^c) = 0, we find 



t{t)= / X X ( -dtA{x,t) + {uj{t) X x) xV X A{x,t))fe{\x\)d^x 

d C C 

= -— / X X A{x,t)fc{\x\)d^x + X X f{u:{t) X x) X V X A{x,t)jfc{\x\)d^x. 
Jr3 Jrs V / 



(5.8) 



The last term in (R.8|) can be rewritten as 



XX {{uj{t)xx)xV xA{x,t))fe{\x\)d^x = uj{t)x / X X A{x,t)fe{\x\)d^x. (5.9) 

To verify (|5l9|) , first note that £c x ((cj(t) x a;) x V x ^(x,t)) = uj{t) x x(^x ■ V x A{x,t)) (for 
a; ■ (lj X a;) = 0) and pull uj{t) x in front of the integral, next use V x a; = and another standard 
identity from vector analysis to rewrite a;-Vx«4 = a;-Vx^ — .A-Vxa; = V-(^xa;), then 
integrate by parts, use the identity (a; x A{x,t)) ■ V/e(|a;|) = and get 

[ a;(a;- V x^(a;,t))/e(|a;|)d^a; = / (a; x ^(a;, t)) ■ V(a;/e(|a;|)) d^x 

= XX A{x,t)fc{\x\)d^x, 
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as claimed. Defining now the electromagnetic field spin vector of the particle by 



Sf(t) = / X X A{x,t)f^{\x\)d^x , (5.11) 

and its canonical spin vector by s = Sb + and finally recalling that a; x s^, = 0, we conclude 
that ([4.5|) can be recast into the canonical evolution equation for the spin (in JFi^b), 



d 

— s = uj X s . (5.12) 
at 



Remark: It follows directly from l\5.1^ that \s\ is conserved during the evolution. 



5.3 The bare spin / angular velocity relation 



Inserting the explicit integral representation for A.{x,t) into the canonical equation ( p. 12 ), and 



recalling that Sb(t) is given in terms of uj{t) by ( [4.3|) , we see that (|5.12|) becomes a closed, non- 
autonomous, nonlinear first-order vector differential equation for a;(t). However, it is advisable 
to eliminate uj{t) in favor of Sb(t). 

We rewrite (O) as Sb(t) = Ib(|ct;(t)|) ■ a;(t), where 



Ib(|a;|)=/ ML^£^/„^(|^|)d3,;, (5.13) 



II |2 

1 — |a; X a;| 

is the inertia tensor of the bare particle. Clearly, lb acts as a number on uj, viz. lb ■ = XbCJ. 
Performing the angular integrations we are left with 

Th{H) = 2nj-^ f^ \ ( [e + 1) artanh(0 - i)idi . (5.14) 

By hypothesis, < Xb(0) < oo. This implies that the map |ci;| i— >• Xb(|a;|) is strictly positive, 
increasing, and strictly convex for |ci;| G [0, 1/-R). Depending on the choice for /m, the bare spin 
magnitude |sb| may or may not approach a finite limit as /" 1. In any event, it follows 
that for |sb| < we can invert the map a; i-^ Sb = Xb(|a;|)a; to get the Euler angular velocity 
vector UJ uniquely in terms of the bare spin vector Sb? viz. uj = 'W(sb), where 

W(sb) = ^(Jbid)-i(|sb|) for |sb|<4. (5.15) 

l^bl 

Note that the map |sb| ^ (Xb id)^^(|sb|) is bounded, strictly increasing, and concave, hence it 
has its steepest slope when |sb| -^0^. This slope at 0"*" is simply the reciprocal value of the slope 
of the map |a;| i-^ |Lj|Xb(|<-jn at \uj\ 0+, viz. slope of (Xbid)~^ ^ Xb(0) (< oo, for Xb(0) > 0, 
by hypothesis). Finally, if sf, < oo, we extend W continuously different iably to M'^ by setting 

W(sb) = l^bl ^4- (5.16) 

it Sb 
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5.4 Bare spin evolution as fixed point problem 

Substituting H'(sb) for in ( |5.12| ) and integrating ( p.l2D w.r.t. t, supplementing the initial 
datum Sb(0) (automatically compatible with the subluminality requirement jtt'ol-R ^ 1), and 
writing out dependencies on A explicitly, we arrive at the following integral equation for Sb, 



Sb(t)=Sb(0)+ / xx[A^{x)-A{x,t)\f,(\x\)d^x 
+ / W(sb(i)) X / a; x A{x,t)f^(\x\)dc'xdt, 

Jo 



(5.17) 



where A = Awavc + ^source is given by the integral representations (|5.7|) and (|5l6|), and where 



uj{t) = W(sb(t)) in (|0|) , closing the chain. Substituting (^]6|) for ^source in ( ^.17|) and rear- 
ranging some integrations gives the explicit fixed point equation for Sb, 

Shit) = Sb(0) + XX (^„avc(a;,0) -^„avc(a;,t))/c(|a;|)d^x 

t 

"/ (>^K(t))-^o)^(t-t)dt 

+ / W(sb(t)) X / x xXavo(ic,t)/o(|ic|)d3xdt (5.18) 
Jo Jr^ 

-Wo X / W(sb(t)) / i^(t')dt'dt 

+ / w(sb(t)) X /" y\;{s^{t'))K{i-t')dt'di, 

Jo Jo 
where K is the electron's retarded self-interaction kernel, 



^(^) = o ] ^U\x\)U\ymt-\x-y\)d'xd'y. (5.19) 

3 Jr3 Ju3 \x - y\ 

Notice that K G L°°(M), and that supp(-ft') C [0,2i?]. By the S0{3) invariance of fe we can 
carry out the angular integrations in ( |5.19|) , obtaining a double integral, 

K{t) = — / Q{t-\r-s\)Q{r + s-t){r'^ + s^ -t^)rsfe{r)f^{s)drds. (5.20) 
3 Jo Jo 

5.5 Lipschitz estimates 

Lemma 1: The map -.M? ^ M? is Lipschitz continuous for the standard Euclidean norm, 
with Lipschitz constant 1/Xb(0). 

Proof. The sole action of W is to scale any input vector u by the factor |'W('u)|/|n|, with 
|>V(n)| = (Jbid)-i(|u|) for \u\ < sf,, and |>V(n)| = 1/R for \u\ ^ sf,- The map \u\ ^ 
(Xb id)~^(|w|) is increasing and concave, vanishing with finite slope 1/Xb(0) for Inl —>■ 0+, and 
saturating for \u\ to (Xb id)~-'^(|n|) 1/-R, with vanishing slope. Thus, |>V(n)|/|'u| is 

monotonic decreasing and bounded above by lim|.„|^o+(|l^('^*)|/|'*^|) = l/2^b(0). Hence, the map 
Xb(0)'W shrinks any input vector it by a factor which is the smaller the longer u is, but leaving 
its direction unchanged. It now follows right away that |Xb(0)>V('Ui) — Xb(0)H'('a2)| ^ I**! — 
for any two vectors Ui and U2. QED 
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Lemma 2: The two-point map : x — > defined by {u, v) y-^ >V(it) x W{v) is 

Lipschitz continuous with Lipschitz constant (Xb(0)-R)~-'^. 

Proof. By a simple identity, followed by the triangle inequality, followed by the upper bound 
\W\ and by Lemma 1, we find 



X W{v^) - W(n2) X W{v2] 



W{ui) X (W{vi) -W{v2)') + (W(ni) ->V(li2)) XW{V2) 



QED 



^ \W{u{} \ \W{V^) -W{V2)\ + \W{ui) - W(W2)| |W(^;2)| 
< ( \W{V,) - W{V2) I + - W{U2) I ) 

^ {RXy,{0)y\\u, - U2\ + \v^ - V2\) 

Writing ( |5.18| ) as Sb = ^(sb) defines a map ^ in the space Lj^(]R+, R^) of M^-valued functions 
u on M"*", equipped with the weighted norm H'uHj^^ = exp{—\t)\u(t)\dt, A > 0. Since by 
assumption the integral of the wave fields over the particle support is bounded for all t, there 
exist two constants Ci and C2, determined by the initial data alone, such that < Ci + 6*2^ 

for any u E L\{R+,R^). Hence, ^ maps Lj^(M+,R3) into some ball {\\u\\^^^ ^ C} C L^(R+,R3), 
where C is determined by the initial data. This also implies that ||.Sb||i a is well defined for any 
solution of ( |5.18| ). 

Proposition 1: The map u J^{u) is \\ . \\^ ^-Lipschitz continuous with Lipschitz constant 



1 



AXb(O) 



' wave 1 1 QQ 



+ 1 



\K\\ 



- II^Im I , 



where ||swavc|L = supt \J^3X x Awaye{x,t)fe{\x\)d^x\ < 00; \\K\ 



(5.21) 

suPtg[o,2i?] < 00, and 



|ir||, = /o^"|ir(t)|dt<oo. 
Proof. By definition of J^, 



J^{u)it) -J^{v){t)= -j^ (W(n(t)) -W{v{t))^K{t-i)dt 

+ l"(w{u{i))-W{v{i))')x I xxA^.ye{x,i)f,{\x\)d''xdt 



'0 

-^0 X 



W{u{i)) -W{v{i))'^ j_ K{t')dt'di 



(W(n(t)) X W{u{t')) - W{v{i)) X y\;{v{t'))^K{i -t')dt' dl 
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Subadditivity of the norm gives 



||:F(n)-^(t;)||,_, 



W{u{i)) - W{v{i))K{t - i)di 



1,A 



+ 



Wo X ^ {W{u{i)) - W{v{i))^J^ K{t')dt'di 
j^j^ (W(u(t)) xW{u{t')) - W{v{t)) xW{v{t'))^K{t-t')dt' di 



1,A 



1,A 



We now estimate one by one the terms on the right-hand side. For the first term we find 
/ (W(u(t)) - W{v{t))K{t - i)di 

Jo ^ ^ 1,A 

|VW(w(t)) - y^{y{t)) I \K{t - i)\dldt 

r iw(w(f)) - I didt 

= II^ILA-W e~''\W{u{t))-W{v{t))\dt 
Jo 

POO 

^ l|i^lloo(A^b(0))-7 e-'' \uit) - v{t)\ dt 

Jo 



^ \\Kl 



ll^lloo (A:^b(0)) 



\u — v\ 



1,A ' 



where in the third step we used integration by parts together with |W| < 1/R and with te 
for t = and t — > 00. The last estimate then is Lemma 1. Similarly, for the second term we 



J^(w{u{i))-W{v{i))) 

poo nt 

/ e-^*/ W{u{t)) -W{v{t)) 
Jo Jo 

poo pt 

/ e-^w w{u{t)) -y\;{v{t)) 

Jo Jo 

W[u{t)) - W{v{t)) 



^ lis 



1,\ 



wave 1 1 00 



Is II X-1 
I •'wave I loo ^ 



X / XX A^^-^e{x,t)f^{\x\)d^xdi 

didt 
dt 



dtdt 



wave 1 1 00 
Swavelloo 



poo 

(A2-b(0))"W e-^^\u{t) - v{t)\dt 
{Xly,{0))-'\\u-v\\,^,. 
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Proceeding analogously for the third term, we find 



X 

< 1^0 



I-2R 







W{u{i)) -W{v{i)) 



2R 



1,A 



^ I'^ol ll^lli 
= \(jL>o\ \\K\\-^ A 



W{u{i)) -W{v{i)) 



K{t')\dt' didt 
didt 



e-^' \W{u{t)) -W{v{t)) \ dt 



e-^' \u{t) - v{t)\ dt 



|a;o| ll^lli (AXb(O)) \\u 



v\ 



1,A • 



For the fourth term we need Lemma 2, otherwise we proceed along the same lines to find 
^ ^ (W{u{i)) X W{u{t')) -W{v{i)) X W{v{t'))^K{i-t')dt'di 

f f \W{u{t)) X W{u{t')) - W{v{t)) X W{v{t'))\\K{i ~ t')\dt' didt 
Jo Jo ^ 

[ \W{u{t)) xW{u{t')) -W{v{t)) xW{v{t'))\\K{t-t')\dt'dt 

° oo t 

^ (Ai?Jb(0))"'^"-^*^ (^\u{t) - v{t)\ + \u{t') - v{t')\^ \K{t-t')\dt'dt 

= (Ai?Xb(0))"Y""^*( - v{t)\ j' \K{t')\ dt'+ y* \u{t') - v{t!)\ \K{t - dt')dt 

^ {\RUO))-'^\\K\\, \\u - v\\,^, + \-' \\K\\^ \\u - v\\,^,y 

Adding all estimates together and finally noting that |a;o| -R ^ 1, we find that 

\\T{u) -J^iv)\\^^^ ^ L\\u-v\\-^^^ 
with L given in (|5^ ). QED 

5.6 Global well-posedness 



The existence of a unique || . ||^ ^^-strong forward solution t Sb(^); t ^ 0, of ( 5.12 ) now follows 
right away from the || . ||-,^ ^-Lipschitz continuity of ^ and the fact that ^ maps Lj^(M'^,M^) 
into some ball || . ^ C*, with C determined by the initial data. Moreover, we can exchange 
t ^ —t and the conclusions holds for the backward evolution as well. Furthermore, for any 
permissible incoming data (not necessarily scattering data) .4wavc(*,0) we can find a A* such 
that L < 1 for all A > A*. We summarize these findings in the following theorem. 
Theorem 1: There exists a unique \\ . \\^ -f^- strong solution t i— > Sb(t) of l \5.12{ ) globally in 
t e M. Furthermore, for all A > A* the map ^ is a \\ . \\^ ^-contraction mapping, and in these 
norms the simple iteration 



'b /' 



(5.22) 



starting with initial datum s^^ = Sb(0), converges || . \\-^^ strongly to the solution t t—>- S\^{t). 
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Global well-posedness in || . ||-,^^ can now be bootstrapped to higher regularity for Sb(t), 
e.g. regularity if the Cauchy data for ^wave and the densities /m and /e are of class Cq(A), 
where A C M'^ is compact. The regularity of A.{x,t) follows accordingly. Unfortunately, a 
detailed discussion of higher regularity is beyond the scope of this letter and has to be deferred 
to some later work. However, note that analyticity of t h- *■ Sb(t) cannot hold, first because of the 
compactly supported ^wave, /m and /e, and furthermore, because it takes only a finite amount 
of energy to spin up the particle so that its equatorial velocity reaches the speed of light when 
/m(| .\)eL'^ (with compact support in Br C R^); for Mb(t) = M\,{\u){t)\), where (by the 
S0{?>) invariance of /m) 

M^{\uj\) = A-K—- / /m ^ artanh(Oede, (5.23) 



and we see that |A^b(|<-<-'|)| < C as |a;| \/R whenever /m G . 

5.7 Conservation laws 
Proposition 2: The following quantities are conserved during the evolution: 

pd^x {charge), (5-24) 



(5.25) 



^ = ^1 {\S\^ + \B\'')d\ + M^{\uj\) {energy), 



C = — [ X X {S X B) d^x + Sb {angular momentum), (5.26) 

a = |sb + Sf| {canonical spin magnitude) . (5.27) 

Proof. We basically follow 0] where the conservation laws for the semi-relativistic theory are 
discussed. 

As for charge conservation, by way of construction LED honors the continuity equation 

dtp{x,t) + V ■j{x,t) = 0, (5.28) 

where p is the electric charge density and j the vector of the electric current density, and this 
fact does not change by simply imposing the condition that the world-line be straight. Indeed, 
one directly verifies that for our j{x,t) = fe{\x\) u!{t) x x we have V ■ {uj{t) x xf^{\x\)) = 0, 
and of course p{x,t) = /e(|a^|) independent of t, i.e. dtp{x,t) = 0. Hence, charge is conserved. 
As for the energy conservation, taking the time derivative of the field energy gives us 



dt V87r 



/ {\S{x,t)\'+\B{x,t)\')d'x) = - [ S{x,t)-j{x,t)d?x, 



(5.29) 



here with j{x,t) = fe{\x\) uj{t) x x. On the other hand, by direct calculation with (|4.7|) and 
one readily verifies that 

^M^{\iv\) = iv^s^. (5.30) 
at at 
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Next, taking the Euchdean inner product with u; on both sides of the canonical evolution equation 
for the total spin, ( ^.121) , we see that 



■ — Sb = -'^ ■ 3-Sf- (5.31) 
at at 

Recalling now the definition of the electromagnetic field spin, ( [5. Ill ), then using the cyclicity of 
a; ■ (a; X dtA)., noting next that ~dtA = S + V^coui and that [oj x x) ■ 4>co\i\{\x\) = 0, and at 
last recalling that fe{\x\) u!{t) x x = j{x,t) we find 



d f 
-oj{t) ■ —st{t) =-Lj{t)- / xxdtA{x,t)f,{\x\)d^x 
dt Jus 



[ {iv{t)xx)-£{x,t)U\x\)d^x= [ £ix,t)-j{x,t)d^x. 



(5.32) 



Hence, energy conservation is proved. 

As for the angular momentum conservation, taking the time derivative of the field angular 
momentum gives the well-known formula H] 



d_^ 1 
dt 



(— [ XX (£(x, t) X B(x, t)) d^x] = - [ XX (p(x, t)£(x, t) + jix, t) x Six, t)) d^x. 

(5.33) 

Inserting our expressions p{x,t) = /e(|a:;|) and j{x,t) = /e(|a;|) uj{t) x x, we see that 

/ X X (p{x, t)£{x, t) + j{x, t) X B{x, t)) d?x = — Sb(t), (5.34) 
Jr3 V /at 

and conservation of angular momentum is proven. 

Finally, we already remarked that ( ^.12| ) implies at once that |s| is conserved. QED 

The proof that the total energy is conserved has the following spin-off. 
Corollary 1: The constraint equation ( U-iA) is automatically satisfied by any solution of 
gyroscopic LED. 

As for the total linear momentum, 

V = ^l £ xBd^x + Ne-uj, (5.35) 

we remark that conservation of ( p. 351 ) is equivalent to ( [4.12[ ). However, our assumption of a 
straight particle world-line is generally not compatible with ([4.12|) , unless special symmetries 
prevail. An example is discussed in the next section. 



6 Rotation-reflection symmetric scattering 

Our Theorem 1 reduces the global existence and uniqueness problem for proper LED with a 
straight particle world-line to finding the class of non-stationary initial conditions for which 
momentum conservation holds with a non-moving particle. Such a class of initial conditions 
is given by the rotation-reflection symmetric field decorations of spacetime, with the particle's 
axis of rotation necessarily identical to the axis of symmetry a, linearly superimposed on which 
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is a compactly supported, non-symmetric electromagnetic radiation field that is never going to 
interact with the particle. Since a non-interacting radiation field is evidently rather uninteresting, 
we confine our discussion to the rotation-refiection symmetric evolutions. 

More precisely, let {(, 6, z) denote cylindrical coordinates of cc, with origin in the particle 
center, axis unit vector a, z = x ■ a, 6 the polar angle of x about a, and ( = \x — za\. The 
axis a is fixed during the evolution, and oj = ua, so that u is (assumed, and below verified 
to be) the only remaining dynamical degree of freedom of the particle. Aside from the non- 
dynamical and spherically symmetric Coulomb field ( ^.5|) , the remaining electromagnetic field is 
now determined by a vector potential of the form A{x, t) = ip{C: t)^^: satisfying the refiection 
symmetry ilj{(,z,t) = il){C,.,—z,t), and obviously rotation invariant. The inhomogeneous wave 
equation for A, (|5.3|) , reduces to the inhomogeneous, scalar, generalized wave equation 



[du - + C^'^c - 5..)^(C,^,t) = 47ra;(t)CVc(V'e+^), (6.1) 

with accordingly simplified scalar solution formulas for ip. An elementary calculation with S = 
— V0coui — dfijj'VO and B = 'Vijj x V6' then shows that the torque J x x {S x B)fed^x oc a, 
establishing the consistency at the level of the gyroscopic problem, indeed. 

6.1 Momentum balance 



We already saw that the time component ( |4.11| ) of the covariant world-line equation is automat- 



ically satisfied, see section 5. We now show that for rotation-refiection symmetric solutions to 
the gyroscopic problem the space-part of the world-line constraint equation ( [4.12| ) is satisfied. 



too. Since the fulfillment of ( [4.12| ) is equivalent to the conservation of linear momentum ( ^.35| ), 
it suffices to show that ( [5. 35] ) is a constant vector for all time. 

By direct computation with S = — V0coui — dtip'Vd and B = x ^6 one verifies that 

/ £{x,t) X B{x,t) fMx\)d'x = - [ C-%^{C,z,t)Vij{C,z,t)d^x = (6.2) 

for our rotation-refiection symmetric fields. As for the spin-orbit coupling term, another direct 
calculation yields that rotation-refiection symmetry implies 

1 p 

N,{t)-iv{t) = -uj{t)— x^{C,z,t)U\x\)d^x = 0, (6.3) 



and the satisfaction of the world-line constraint equation ( [4.12[ ) follows. 



6.2 Exponential convergence to the soliton state 

In m we proved that the conservation of cr = |sb + Sf| together with the invertibility of the 
map a; I— >■ s in stationary situations implies that any scattering process connects two boosted 
stationary particle states with identical values for the renormalized mass and the magnitudes 
of spin and magnetic moment. In short: the Lorentz electron scatters like a soliton. We now 
complement this result by proving that rotation-refiection symmetric scattering does occur, and 
that the soliton state is approached exponentially fast. For our proof we need to assume that 
the ratio of electrostatic to bare rest mass is small. 
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Proposition 3: Assume that the electromagnetic potential data are rotation-reflection sym- 
metric in the sense explained above, and of class C^. Assume furthermore that ipwa.ve{C, z,0) has 
compact support a finite distance away from supp(/e). Finally, assume that 

ll/^lli < Jb(0). (6.4) 

Then, as t ^ oo, the bare spin Sb(t) converges exponentially fast to a stationary vector, Sb(^) 
s^, and = Sb(0). 

Proof. Clearly, since oc a for all t, all terms c^o x W(sb) and W(sb(^)) x W(sb(t)) vanish. 
Also, by direct calculation one verifies that J^^x x Awa.ve{x, t)fe{\x\) d^x oc a for all t, so that its 
cross product with W vanishes as well for all t. Furthermore, by hypothesis, the initial wave data 
don't overlap with the support of the particle, hence J^^x x Awstve{x, 0)/e(|a;|) d^x = 0. Finally, 
by the wave propagation, there exists a T ^ 2R such that supp [Avawd^, t)) flsupp (/c(|a;|)) = 
for all t > T. Then, for t > T, we have 

Sb{t)+ [ W{sdt))K{t-i)di= s{0), for t>T (6.5) 

Jt-2R 

where s(0) = Sb(0) + kojq, with k = ^_ft'(t)dt. Notice that ( |6.5|) is effectively a scalar equation 
because all vectors are oc a. We now define s'^ as the - unique - solution of 

s^ + kW{s^)=s{0). (6.6) 

Clearly, since s(0) = Sb(0) + kojq and >V(sb(0)) = u^o, (§]§) is solved by = Sb(0), and by 
uniqueness this is the only solution. We next rewrite (|6.5|) as 



sdt)-s^ = -l" (w{sdi))-W{s^))K{t-i)di for t>T 



(6.7) 



and estimate 



isb(t)-ci ^ r iw(sb(t))-w(sc°)iiir(t-t)idt 

Jt-2R 

^ (Xb(0))"' r \sdi)-s^\\K{t-i)\di (^-^^ 

Jt-2R 

^ (^b(0))"VaXte[t_2i?,i] \Mi) - ^h] 

where we used the Lipschitz continuity of VV (Lemma 1) and the continuity of t i— > S]^{t). Now 
assume that t G [n2R, (n + l)2i?], with n big enough so that n2R > T. By ( |6.8| ) and the inclusion 
[t - 2R, t] C [{n -l)2R,{n + 1)2R] we have that 

max |sb(t) - ^ IlirlL (Jb(0))'^ max |sb(t) - sS"! (6.9) 

te[n2R,{n+l)2R] ^ ^ ^ te[{n-l)2R,{n+l)2R] ' ^ ' " ' 

By our smallness condition ( |6.4|) we conclude that maxjg[(„_i)2R,(„+i)2_R] \sh{t) — cannot be 
attained in [n2R, {n + 1)2/2], hence it is attained in [{n — l)2R,n2R). By induction from one 
interval of length 2R to the next one we now get \si,{nT) — ^ C exp(— nF), i.e. exponential 
convergence with rate F = In (Xb(0)/ \\K\\^ ). QED 



The exponentially fast convergence Sb(t) implies for all rotation-reflection symmetric 

initial conditions of the type discussed that the field-particle system in fact converges expo- 
nentially fast on families of nested compact sets to a stationary particle-field bound state, the 
soliton state, while a departing field of electromagnetic radiation escapes to spatial infinity. Put 
differently, our class of rotation-reflection states consists of scattering states, with the exception 
of the stationary bound state itself. For late times the evolution of the electromagnetic field thus 
satisfies the scattering formulas (* means complex conjugate) 

G{x,t) Q::^{x) + e-^'^^gZ{^), (6.10) 

and 

gZ*{x) + e^'^-gil^*{x), (6.11) 
where the soliton fields and ^""j* coincide in this rotation-reflection symmetric setting. 



7 Open problems 

It is instructive to have some explicit numbers. As in [|I|], consider the example where /e 
and /m are given by the uniform surface measure on a sphere of radius R, i.e. /e(|a?|) = 
—e{A7cR'^)^^ 5{\x\ — R), and /m(|a;|) = mb(47r_R^)~^ (5(|a:;| — R), with the strictly positive 
bare rest mass of the Lorentz electron. This gives Xb(0) = (2/3)mb-R^, and 

m = e'^-{l-^^)e{t)e{2R-t), (7.1) 

so that \\K\\-^ = e^i?2(2-\/2 — l)/9. Our smallness condition \\K\\-^ < Xb(0) then becomes 



m^R 2V2-I' 

Roughly speaking, the particle's electrostatic Coulomb energy must be less than the bare rest 
mass. (This conclusion holds with minor numerical differences also when /m is uniform volume 
measure in Br.) The interesting question now is whether deviations from the soliton state decay 
exponentially fast also when the smallness condition ( [7.2|) is violated, especially since one is 
interested in a renormalization flow limit —* 0^ where R 1.5Rc (with R^ the electron's 
Compton length) |l|. Conceivably some long-lived resonances may emerge and render a more 
complicated picture. Nonlinear resonances have been studied rigorously in the simpler semi- 
relativistic model of a particle interacting with a scalar wave field |Q; see also |jl3| for certain 
nonlinear wave equations. A corresponding study for gyroscopic LED is in its infancy. 

For general non-rotation-refiection symmetric initial data we proved global existence and 
uniqueness of gyroscopic solutions (which typically do not satisfy the world-line equations of 
LED), but we do not yet know that on families of nested compact sets the field-particle system 
converges to a stationary state. All we can show is that Sb(^) converges to some as t ^ 00 

whenever the iterated integral /q /q l^(^b(^)) x {shit')) K {i — t')dt' di has a limit in as 
t — >■ 00, but we have nothing to say about exponentially fast convergence, then. In case of 
a scattering scenario, i.e. with convergence to a soliton, the fields and are generally 
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not identical; however, they differ by at most a space rotation as a consequence of the sohton 
dynamics. The exphcit characterization of the scattering operator from the "in" states to the 
"out" states has yet to be worked out. 

Eventually we would like to be able to establish control over the problem of many-body scat- 



tering. While well developed in quantum theory g, ^, |TT], 0], very little is known rigorously 
for truly relativistic LED. Interestingly enough, the solution to this problem requires the con- 
struction of a self-consistent nontrivial foliation of space-time, injecting a technical element from 
general relativity into the analysis. 
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